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Abstract 

Oliinychenko, Bugaev and Sorin [arXiv: 1204.0103 [hep-ph]] considered the role of conservation laws 
in discussing possible weaknesses of thermal models which are utilized in describing the hadron multi- 
plicities measured in central nucleus-nucleus collisions. They argued to analyse the criteria for chemical 
freeze-out and to conclude that none of them were robust. Based on this, they suggested a new chem- 
ical freeze-out criterion. They assigned to the entropy per hadron the ad hoc value 7.18 and supposed 
to remain unchanged over the whole range of the baryo-chemical potentials. Due to unawareness 
of recent literature, the constant entropy per hadron has been discussed in Ref. [Fizika B18 (2009) 
141-150, Europhys.Lett. 75 (2006) 420]. Furthermore, it has been shown that the constant entropy 
per hadron is equivalent to constant entropy normalized to cubic temperature, an earlier criterion for 
the chemical freeze-out introduced in Ref. [Europhys.Lett. 75 (2006) 420, Nucl.Phys.A764 (2006) 
387-392]. In this comment, we list out the ignored literature, compare between the entropy-number 
density ratio and two criteria of averaged energy per averaged particle number and constant entropy 
per cubic temperature. All these criteria are confronted to the experimental results. The physics 
of constant entropy per number density is elaborated. It is concluded that this ratio can't remain 
constant, especially at large chemical potential related to AGS and SIS energies. 
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I. INTRODUCTION 



In the preprint 1|, Oliinychenko, Bugaev and Sorin have considered the role of conservation 
laws, the values of hard core radii along with the effects of the Lorentz contraction of hadron 
eigen volumes in discussing the weaknesses of thermal models which are utilized in describing 
the hadron multiplicities measured in the central nucleus-nucleus collisions. Regardless the 
unawareness of earlier literature, the authors concluded that none of the criteria for the chem- 
ical freeze-out is robust. In doing this, they entirely disregarded the experimental results in 
baryo-chemical potentials /if, and their corresponding temperatures T. A systematic analysis 
of the four criteria describing the chemical freeze-out is introduced in Furthermore, a 

comparison between these four criteria is elaborated in {2 -4 1. 

Starting from phenomenological observations at SIS energy, it was found that the averaged 
energy per averaged particle (e) /(n)«l GeV 



5[, where Boltzmann approximations are applied 



in calculating (e) / (n) , this constant ratio is assumed to describe the whole T — fif, diagram. 
For completeness, we mention that the authors assumed that the pions and rho-mesons get 



dominant, at high T and small fih. The second criterion assumes that total baryon number 
density (n^) + (ni) ~ 0.12 fm -3 |6|. In framework of percolation theory, the authors of Ref. 



3[, the last two criteria seem to 



have suggested a third criterion. As shown in Fig. 2 of 
give almost identical results. All of them are stemming from phenomenological observation. A 

|3j. Accordingly, the 



fourth criterion based on lattice QCD simulations was introduced in Ref. {2], 
entropy normalized to cubic temperature is assumed to remain constant over the whole range 
of baryo-chemical potentials, which is related to the nucleus-nucleus center-of-mass energies 
\/snn Q- An extensive comparison between constant (e)/(n) and constant s/T 3 is given in 

The thermodynamic quantities deriving the chemical freeze-out. In framework of hadron 
resonance gas are deduced {2, 3|. Explicit expressions for s/n at vanishing and finite temper- 



ature are introduced [2j, [8| . The motivation of suggesting constant normalized entropy is the 
comparison to the lattice QCD simulations with two and three flavors. We simply found the 
s/T 3 = 5 for two flavors and s/T 3 = 7 for three flavors. Furthermore, we confront the hadron 
resonance gas results to the experimental estimation for the freeze-out parameters, T and \x\,. 



2 



II. THE HADRON RESONANCE GAS MODEL 



The hadron resonances treated as a free gas |s-ll3| are conjectured to add to the thermo- 
dynamic pressure in the hadronic phase (below T c ). This statement is valid for free as well as 
strong interactions between the resonances themselves. It has been shown that the thermody- 
namics of strongly interacting system can also be approximated to an ideal gas composed of 
hadron resonances with masses < 2 GeV jl2, [l4]. Such a mass cut-off is implemented to avoid 



the Hagedorn singularity 15|. Therefore, the confined phase of QCD, the hadronic phase, is 
modelled as a non-interacting gas of resonances. The grand canonical partition function reads 

Z(T,V) = Tr[exp- H / T ], (1) 

where H is the Hamiltonian of the system and T is the temperature. The Hamiltonian is 
given by the sum of the kinetic energies of relativistic Fermi and Bose particles. The main 
motivation of using this Hamiltonian is that it contains all relevant degrees of freedom of 
confined and strongly interacting matter. Obviously, it can be characterized by various - but 
a complete - set of microscopic states and therefore the physical properties of the quantum 
systems turn to be accessible in approximation of non-correlated free hadron resonances. Each 
of them is conjectured to add to the overall thermodynamic pressure of the strongly interacting 
hadronic matter. It includes implicitly the interactions that result in resonance formation. 
In addition, it has been shown that this model can submit a quite satisfactory description of 
particle production in heavy-ion collisions (q- 13] . With the above assumptions the dynamics 
the partition function can be calculated exactly and be expressed as a sum over single-particle 
partition functions Z\ of all hadrons and their resonances. 



In Z(T,^,V) = ^ln^(T,y) = ^±g^ °A; 2 ^ln(l±exp[^^^ 



(2) 



where €i(k) = (k 2 + m 2 ) 1//2 is the i— th particle dispersion relation, is spin-isospin degeneracy 
factor and ± stands for bosons and fermions, respectively. 

The switching between hadron and quark chemistry is given by the relations between the 
hadronic chemical potentials and the quark constituents; /i, = 3nt,fi q + n s fis, where nb(n s ) 
being baryon (strange) quantum number. The chemical potential assigned to the light quarks 
is \i q = (fi u + fid)/ '2 and the one assigned to strange quark reads Us = \i q — \i s - The strangeness 
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chemical potential fis is calculated as a function of T and fii under the assumption that the 
overall strange quantum number has to remain conserved in heavy- ion collisions 12]. 

The HRG calculations assume quantum statistics and an overall strangeness conservation. 
With this regard, the strangeness chemical potential /15 is calculated at each value of T and fib 
assuring that the number of strange particles should be the same as that of the anti-strange 
particles. It is worthwhile to mention that no statistical fitting has been applied in determining 
all thermodynamic quantities, including entropy and number density derived from Eq. (J2J). 

III. PHYSICS OF CONSTANT ENTROPY PER NUMBER DENSITY 

From the entropy and equilibrium, the Gibbs condition simply leads to 



the rhs is positive as long as /i& < p/n + e/n, where the thermodynamic quantities, p, e and n 
are supposed to be calculates at the T — fib diagram of the chemical freeze-out. Fig. [1] shows the 
experimental estimation for the freeze-out parameters T and fib- It is obvious that increasing 
fib leads to decreasing T and therefore all values of the thermodynamic quantities decrease as 
well. Cleymans et al. [s| suggested an empirical T — fib relation 



where a, b and c are fitting parameters. In light of this discussion, the value given to s/n can't 
remain unchanged with increasing fib- Left panel of Fig. [2] presents the values of the three 
criteria (e)/(n), (rib) + and s/n calculated in HRG, section HT1 at s/T 3 = 7. It is obvious 
that a„ four criteria S ee m to remain constant, espectaUy at h lg h At low energies, the 

value assigned to s/n [1[ is larger than the actual one, the value resulted from order conditions. 
The reason is illustrated in the right panel. At different values for fib, the thermal evolution 
of s/n is presented. It is obvious that s/n never reaches 7.18 at fib > 500 MeV. It is essential 
to bear in mind that the value 7.18 has almost no physical interpretation. It is just an ad 
hoc value. This makes it inapplicable at Alternating Gradient Synchrotron (AGS) and Schwer- 
Ionen- Synchrotron (SIS) energies. Almost same kind of restriction would be valid for e/n. 
According to Eq. ([3]), 




(3) 



T = a-bfi 2 b -cfi 4 b , 



(4) 



The physics of constant s/T 3 has been discussed in Ref. j^, 3]. It combines the three 
thermodynamic quantities, p/T A , e/T 4 and n/T 3 

rp 3 + j, 4 H-b rp 3 ■ {O) 

At chemical equilibrium, the particle production at freeze-out is conjuncted to fully fulfil the 
laws of thermodynamics, as Eq. ([3]). The hadronic abundances observed in the final state of 
heavy- ion collisions are settled when s/T 3 drops to 7 i.e., the degrees of freedom drop to 7n 2 /4. 
Meanwhile the changing in the particle number with the changing in the collision energy is given 
by Hb, the energy that produces no additional work, i.e. the stage of vanishing free energy, gives 
the entropy at the chemical equilibrium. At the chemical freeze-out, the equilibrium entropy 
represents the amount of energy that cant be used to produce additional work. In this context, 
the entropy is defined as the degree of sharing and spreading the energy inside the system that 
is in chemical equilibrium 



IV. CONSTANT ENTROPY PER NUMBER IN LATTICE QCD SIMULATIONS 
AND HEAVY-ION COLLISIONS 

Once again, related literature on lattice QCD simulations is not cited in For example, 
Borsanyi et al. [16| studied the trajectories of constant s/n, where s = S/V and n = N/V, 
on the phase diagram and thermodynamic observables along these isentropic lines. This was 
not the only work devoted to such line of constant physics 17|. In Stefan-Boltzmann limit, the 
ratio s/n is assumed to remain unchanged with increasing ^ (Appendix A of 16j). In doing 
this, lowest order in perturbation theory is assumed, where strangeness chemical potential us 
likely vanishes. For fib/T, a limiting behavior for the isentropic lines on the phase diagram is 



obtained. The ratio s/n has been measured at various ^snn [IS 



It is concluded that in limits 



of low temperatures, increasing the chemical potential results in an overestimation for the ratio 
s/n even beyond the applicability region of the Taylor-expansion method, which is applied in 
lattice QCD simulations at finite chemical potential. Two remarks are now in order. First, 
the values of s/n seem to depend on the chemical potential /i& or y / 5jvjv. This is confirmed in 
different experiments [18| and lattice gauge theory [16j. Second, the ratio s/n as calculated in 
the lattice QCD simulations jl^] is suggested to characterized the QCD phase diagram 12]. 
The QCD phase diagram is likely differs from the freeze-out diagram 0, 3], especially at large 



chemical potential /i& or small y/JJm so that at fixed //& the critical temperature differs from 
the freeze-out temperature. 



V. RESULTS AND CONCLUSIONS 
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Fig. 1: The freeze-out parameters, T and fib, measured in various heavy- ion collisions experiments 
(labelled) are compared with the three criteria, (e)/(n) = 1 GeV (dashed line), s/n = 7.18 (dotted 
line) and s/T 3 = 7 (solid line). 




In Fig. [U the freeze-out parameters, T and fib, measured in various heavy-ion collisions 
experiments are compared with the three criteria, (e)/(^) = 1 GeV (dashed line), sin = 7.18 
(dotted line) and s/T 3 = 7 (solid line). The experimental data are taken from |4| and the 
reference therein. The quality of each criterion is apparent. All conditions are almost equivalent 
at very high energy. The ability of the condition (e) / (n) = 1 GeV at very low energies are not 
as much as that of s/T 3 = 7. As discussed in section HIT] s/n = 7.18 seems to fail to reproduce 
the freeze-out parameters at fib > 500 MeV. To illustrate the reason for this observation, the 
thermal evolution of s/n at very high chemical potential calculated in HRG is presented in 
the right panel of Fig. |2j Details on HRG are elaborated in section [TXJ It is obvious that the 
value assigned to s/n would be achieved at fib > 500 MeV. In other words, it is obvios that the 
behavior of s/n is non-monotonic. 

The left panel of Fig. [2] presents the energy scan for the three criteria, (e)/(n), (jib) + (ni) 
and s/n calculated in HRG, section HT1 at s/T 3 = 7. The calculations in HRG are performed 
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as follows. Starting with a certain the temperature is increased very slowly. At this value 
of hi, and at each increase in T, the strangeness chemical potential us is determined to assure 
strangeness conservation. Having the three values of fib, T and fis, then all thermodynamic 
quantities are calculated. When the ratio s/T 3 reaches the value 7, then the three quantities 
(e)/(n), (rib) + (rij) and s/n are registered. This procedure is repeated over all values of «&. 
We find that s/T 3 = 7 assures s/n = 7.18 and (e)/(n) = 1 GeV at small fi^ (large ^snn)- 
At large fib (small ^snn), the values of s/n gets smaller values, so that the applicability of 
s/n = 7.18 is limited to fib < 500 MeV. In conclusion, the robustness of s/n = 7.18 is very 
much limited in comparison to the four criteria: percolation 3], baryon number jg], energy 
per particle s| and normalized entropy 

SB. 

That s/T 3 is accompanied with constant s/n 
has been introduced in Ref. That authors of |lj argue that s/n = 7.18 is novel likely 

reflects an ignorance of related literature. The four criteria are based on physical 

observation either phenomenological and/or theoretical. The authors of [l| are suggesting an 
ad hoc value for the ratio s/n. It is inapplicable at AGS and SIS energies. Its relation to 
s/T 3 is apparently overseen. The same is valid for the comparison with other criteria (some 
of them are ignored, completely) and ignorance of the experimental measurements. The ad 
hoc value assigned to s/n is obviously much robuster than any other criterion. Unawareness 
of literature and underestimating or even ignoring previous work are violation of rules of the 
scientific research. 
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Fig. 2: Left panel: the energy scan for the three criteria, (e)/(n), (rib) + (ni) and s/n are calculated in 
HRG, section |TI] at s/T 3 = 7. Right panel: the thermal evolution of s/n at ut> = 500, 630 and 750 MeV 
is presented. The horizontal dashed line indicates s/n = 7.18. The singularities at low temperatures 
are stemming from almost vanishing number density (compare with left panel). 
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